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FOREWORD 


The  nde  that  the  wavelet  transfonnation  mi^t  jday  in  the  identification  of  helicopters  from 
their  acoustic  signals  is  explored  in  diis  document  The  utility  of  using  a  feature  extraction  system 
employing  bqdi  conventional  Fourier  transform  techniques  in  combination  with  the  wavelet 
transform  is  ^scussed.  Depending  on  the  nature  of  fee  mission,  such  a  system  could  mfeer 
improve  fee  helicopter  identificatitm  process  by  suppressing  transient  imise  prior  to  feature 
extraction  using  conventional  techniques  or  piovi^  the  <^>portunity  to  use  transient  signatures  to 
identify  a  qsecU^  helicopter  in  a  groiq)  consisting  (tfhomo^ous  model  types. 

This  work  was  done  in  support  of  fee  Advanced  Processors  for  Weapons  Sensor  Fusitm 
program  of  fee  U.  S.  Marine  and  has  been  conducted  in  fee  Advanced  Computation 
Technology  Group  of  the  Systems  Research  and  Technology  Department 

This  rq)ort  has  been  reviewed  by  Dr.  Richard  A.  Lorey,  Head,  Advanced  Computation 
Technology  Group. 


Af^ffovedby: 

D.B.  COLBY,  Head 

Systems  Research  and  Technology  Dq)aiiment 
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ABSTRACT 


The  ability  to  classify  hdiooplers  baaed  oo  dwir  acousdc  signature  has  ai^licatioas  to  bodi 
detecdon  and  traddng.  The  cmrait  qtproach  to  diis  proUem  uses  standard  signal-processing 
techniques  to  extract  features  baaed  on  ^  fimdamental  hannooics  of  die  helicopters*  blades.  This 
document  examiiies  the  role  that  the  wavelec  ttansfonnadon  nd^  play  in  this  feature  extraction 
process.  It  addresses  ways  to  use  die  wavelet  transfonn  bodi  to  improve  the  performance  of  die 
existing  iqyP*^*P<^>o<^'q>”*^^capd)flides  of  dredassificarion  system 
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INTRODUCTION:  HEUCOPTER  IDENTinCATION  PROBLEM 

The  problem  ci  interest  here  is  die  classification  and  tracking  oi  helicoptm  from  their 
acoustic  signatures.  These  tasks  are  to  take  place  under  field  conditions.  The  timeline  requirement 
for  these  tasks  is  as  close  to  real-time  as  po^le. 

There  are  several  tilings  that  act  to  cooqilicate  tiiis  problem.  First,  tiiere  is  die  existence  of 
transient  batdefield  signals  such  as  gunshots,  explosions,  etc.  Secmid,  there  are  periodic  processes 
tiiat  are  part  al  the  environmoit  These  include  generatin'  noise  and  friendly  cr^  signatures.  The 
final  factor  is  diat  there  are  usually  several  crafts'  signals  acting  in  ctmceit  with  one  another. 
These  frictors  act  either  to  degrade  tte  signal  of  interest  (n  to  introduce  artifacts  inti)  the  signal  that 
copyromise  the  classificatitm  and  tracking  mission. 

Current  tq)proaches  to  this  problem  use  spectral  aiudysis  techniques  to  extract  the  funda¬ 
mental  frequencies  of  the  main  and  tail  rottn  of  the  craft  A  classificatitm  signature  is  formed 
using  the  ratio  of  diese  frequencies.  Recent  work  by  Poston*  et  al.  has  enqiloyed  harmogiam 
analysis  to  extract  these  frequencies. 

This  document  examines  die  use  of  the  wavelet  transformation  (WT)  as  an  alternate 
approach  to  this  feature  extraction  problem.  Feature  extraction  schemes  that  only  use  die  WT  and 
hybrid  wavelet  transformation/fast  Fourier  transfimn  (WT/FFT)  schemes  are  di^ssed  within, 
liie  use  of  the  WT  to  circumvent  some  of  the  above  mentioned  problems  is  also  discussed. 

This  document  begins  with  a  brief  review  of  the  mathematical  underpinnings  of  the  WT, 
including  discussicMis  of  both  cmitinuous  and  discrete  WT.  With  this  background  in  hand,  the 
applicatitm  cf  the  WT  to  signal  processing  will  be  crmsidered  with  specific  attention  to  the  heli- 
c(q;)ter  problem.  These  discussirxis  include  smne  information  about  available  hardware  for  real¬ 
time  inqilementation  of  the  WT.  This  document  concludes  with  recommendations  for  future  work. 


WAVELET  TRANSFORMATION 


As  a  starting  pmnt  for  the  introducticm  to  the  WT,  the  concq>t  of  Fourier  series  expansion 
will  be  reviewed  bricfiy.  In  Fourier  series  expansion,  the  sine  and  cosine  functicms  are  used  at  dif¬ 
ferent  firequencies  as  a  basis  foe  the  qiace  L^(R)  of  square  integrable  functions.  This  basis  is  very 
useful  for  decomposing  functions  of  a  periodic  nature;  i.e.,  those  functions  that  are  localized  in 
frequency  qiace. 

The  basis  functions  used  in  the  WT  ate  dialations  and  translations  of  a  single  function 
known  as  the  mother  wavelet  'P.  The  kth  order  transtation  of  Y  at  a  resolution  of  order  j  is  given 
by 
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(1) 

Given  certain  properties  on  this  mother  wavelet,  the  ^i^’s  form  an  orthonormal  t»asis  of  L^(R). 

A  classical  exanq;>le  of  such  a  modier  wavel^  is  the  Haar  function  yfx) »  1,  if  0  ^  x  ^1/2, 
-1,  if  1/2  ^  ^1,  and  0  otherwise  (see  Figure  1).  The  Haar  function  is  a  good  example  how  sim¬ 
ple  a  mother  wavetetcan  be.  The  wavelets  that  are  used  within  this  document  are  smoother  gener- 
alizatioos  die  Haar  wavelet 


For  the  purposes  of  this  mathematical  develqiment  discussions  will  be  limited  to  the  one- 
dimensitmal  WT.  The  inner  product  of  f<x)  with  g(x)  is  written  as 


{g  («)  ,/(n)  >  =  J  g  iu)fiu)  du 


(2) 


Befme  the  implementadon  details  of  the  discrete  wavelet  transform  (DWT)  are  discussed, 
some  background  cm  the  ctxitinuous  wavelet  transform  is  appropriate.  Given  a  function  f(x),  the 
goal  is  to  eiqiand  the  function  in  the  wavelet  basis.  One  desires  to  write 

/W  =  (3) 
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where  the  expansion  coefficients  are  given  by 

djt=  <f(x),y¥jk)  W 

In  die  DWT  case,  <Hie  starts  with  a  data  set  »  (x^i,  x^2>  •—  where  N  is  a  power  of 
2.  This  discussion  the  DWT  will  be  based  on  the  multiresolution  approach  of  Mallat  In  this 
approach,  one  starts  with  a  blurring  function  9  that  along  with  its  dialations  and  translations  9jk 
spans  the  space  of  multiresolutitm  qipioximations  of  the  original  functk».  This  9.  in  turn, 
induces  a  wavelet  yf  and  basis  i|rjk.  Th^  9jk  and  9jk  act  together  to  span  the  space  of  L^R).  The 
blurring  function  for  the  Haar  wavelet  is  just  the  cluuacteristic  function  <hi  [0,1],  Le.  XOO  -  1  if  x 
e  [0, 1]  and  »  0  otherwise  (see  Figure  2). 


HGURE  2.  BLURRING  FUNCTK^  CXDRRESPONDING  TO 
HAAR  MOTHER  WAVELET 


In  this  manner,  the  original  signal  at  a  resolution  of  2^  is  split  into  a  direct  sum  of  a  multiresolu- 
ticHi  iq[)proximati(Hi  at  a  level  of  2*^  (obtained  using  the  appropriate  4!]k’s)  and  a  detail  signal  (or 
set  of  wavelet  coefficients)  is  obtained  using  the  appropriate  Therefore,  one  symbolically 
writes 


S^  =  S^‘eD^’‘  (5) 

where  is  the  detail  signal  at  a  resolution  of  2'^  In  some  sense,  '  can^be  though  of  as  the 
best  approximation  at  a  resolution  level  of  2'^  of  the  original  signal  and  is  the  set  of  details 
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tfutt  ae  kut  in  goii^  firom  !r  to  S*  .  This  process  can  be  repeated  using  die  multireatdutioo 
approxiination  to  die  signal  a  the  levd  of  2*^  as  the  starting  point  This  leads  to 

=  S®*eD^**  (6) 

which  is  analogous  to  Equation  S.  Finally  we  may  combine  these  equati(»is  together  to  write 

(7) 

In  diis  manner,  die  signal  is  reconstructed  from  its  set  of  wavelet  coefficients  and  lowest  level 
multirescdotion  tqifMoximadon. 

An  inherent  drawback  to  the  Haar  wavelet  is  its  lack  of  smoothness.  This  problem  has 
been  addressed  by  Daubechies^  and  others  by  taking  convolution  powers  of  the  Haar  mother 
wavelet  The  graphs  of  the  Ns2  second'order  convolution  powers  for  the  blurring  function  and 
the  mother  wavelet  are  given  in  Figures  3  and  4. 

By  varying  N,  one  can  obtain  different  blurting  functions  and  wavelet  functions  NT 
The  smoothness  of  these  functitms  increases  as  N  is  allowed  to  increase.  For  a  delightful  introduc¬ 
tion  to  the  creatitm  of  wavelet  bases,  see  the  recent  paper  of  Strichartz.^ 


nOURE  3.  GRAPH  OF  DAUBECHIES  COMPACTLY  SUPPORTED 
BLURRING  FUNCTION  2* 
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-1  0  t  t 


FIGURE  4.  GRAPH  OP  DAUBECHIES  COMPACTLY  SUPPORTED 
MOTHER  WAVELET  FUNCTION  2'*' 


RESULTS 

In  dus  secticm,  some  preliminaiy  classification  results  are  presented.  These  results  were 
olRained  using  field  collected  signals  firom  two  different  helicopter  models.  Ultimately;  the  fielded 
system  must  possess  the  ability  to  distinguish  between  several  different  models,  but  diis  proof  oi 
concept  study  <»ly  addresses  die  two-class  problem. 

The  two  helicopter  signals  are  designated  110  and  117.  The  signals  were  collected  at 
2000  Hz.  The  time  dcnnain  representadtm  for  approximately  2  sec  (4096  samples)  of  signal  data  is 
given  in  Hgure  S. 

The  dyadic  Daubechies  Ns2  wavelet  transform  of  each  of  the  two  signals  is  dven  in  Hg- 
ures  6  and  7.  The  coefficients  are  plotted  bottom  to  top  fixnn  a  resolution  level  of  2*^=  2048  val¬ 
ues  to  a  resolution  level  of  2'^^  »  2  values.  The  decision  to  use  the  N=2  mother-wavelet  for  the 
subsequent  discriminant  analysis  was  made  somewhat  arbitrarily.  Figure  8  shows  the  histogram 
for  the  two  helicq)ters. 
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HGURE  5.  TIME  DOMAIN  PLOT  OF  TWO  SIGNALS 
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FIGURE  6.  WAVELET  TRANSFORM  FOR  HEUCOPTER  110 
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Wtowtol  OaoonipMMon  CotHdanM 


HGURE  7.  WAVELET  TRANSFORM  FOR  HEUCOPTER  117 


•4  M  *.J  *M  *»  1.0 


FIGURE  8.  HISTOGRAM  OF  FEATURES  AT  RESOLUTION 
LEVEL  OF  2*’  FOR  TWO  HELICOPTERS 


Given  this  set  of  coefficients,  one  is  faced  with  the  daunting  task  of  extracting  an  tqipro- 
piiate  set  of  features  from  them.  The  curse  of  dimensionality  and  conunon  sense  prevents  one 
from  conducting  discriminant  analysis  in  With  this  end  in  mind,  the  following  procedure 
was  used. 
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Given  D* « {dV  d2*  ••••  the  set  <tf  wavelet  co^dents  at  level  2^  the  oiergy  E*  at 

this  lesdutioa  level  can  be  estimated  using 

fi*-  X  (») 

i=  1 

The  set  of  12  values  {E^>E^  ....  }  was  then  normalized  by  using  a  standaid  linear  transfonna- 

tion  to  map  the  largest  value  to  1  and  the  smallest  value  to  0.  The  sampling  window  was  then 
shifted  10  units  to  the  right,  and  the  process  was  repeated.  In  this  manner,  591  twelve-dimensional 
data  points  were  generated  for  each  of  the  two  signals. 

These  data  will  now  be  analyzed.  Table  1  contains  the  first  two  feature  vectms  fiv  each 
class.  One  notices  that  the  maximum  energy  concentration  occurs  at  a  different  resolution  level 
for  the  two  classes.  In  particular.  Signal  110  obtains  its  maximum  at  a  resolution  level  2'^, 
while  Signal  1 17’s  maximum  occurs  at  level  2*^.  It  would  be  very  fortuitous  from  a  discriminant 
analysis  perspective  if  diis  trend  persisted  throughout  the  signal  duration.  In  fact,  the  signals  stu¬ 
died  indicate  this  trend  is  consistent  throughout  the  duration  of  the  signal.  The  histogram  of 
Figure  8  illustrates  that  there  is  no  overlap  between  the  histograms  of  the  level  2*^  feature  fw  110 
and  117.  In  fact,  as  with  the  first  two  examples,  the  value  holds  constant  at  1;  i.e.,  the  maximum 
eneigy  for  Signal  1 10  consistently  occurs  at  a  resolution  level  of  2‘^. 

TABLE  1.  FIRST  TWO  FEATURE  VECTORS 
FOR  EACH  CLASS 


Level 

UO-1 

UO-2 

U7-1 

U7.2 

0.000 

0.001 

0.000 

0.002 

1  2-» 

0.004 

0.000 

0.000 

0.000 

0.024 

0.006 

0.012 

0.015 

2^ 

0.253 

0.283 

0.011 

0.011 

2'* 

0.109 

0.071 

0.049 

0.049 

2'^ 

0.137 

0.120 

0.060 

0.057 

2* 

0.344 

0.340 

2’ 

1.000 

0.617 

0.607 

2-* 

0.S96 

0.613 

0.296 

0.271 

2'^ 

0.161 

0.137 

0.118 

0.126 

2*2 

0.106 

0.104 

0.093 

0.085 

2* 

0.111 

0.101 

0.131 

0.129 

It  is  certainly  appropriate  at  this  point  to  elaborate  on  the  significance  of  this  phentmiena 
and  to  try  to  provitte  some  explanation  of  the  und^Iying  mechanism  that  produces  it  Note  that 
perfect  classification  results  can  be  obtained  using  the  firatures  of  Figure  8.  A  P(CC)=1.0  can  be 
obtained  with  a  corresponding  P(FA)sO.O  using  a  simple  threshold  between  the  two  distributions, 
which  cleariy  demonstrates  the  utility  of  this  feature. 
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Given  the  periodic  nature  of  the  signals,  it  is  possible  for  the  energy  to  be  partially  concen¬ 
trated  at  a  particiUar  resolution  level.  Similariy,  the  presence  of  transients  in  the  signal  can  be 
manifested  by  energy  concentration  at  a  certain  resolution  level.  The  phenomena  seen  may  be  a 
ctunbination  of  these  two  effects.  If  the  concentration  is  solely  frequency  based,  then  this  signa¬ 
ture  may  be  useless  in  distinguishing  between  two  of  the  same  models  or  between  any  two  signals 
with  similar  frequency  components.  The  utility  of  the  signature  is  greatly  enhanced  if  it  is  based 
on  transients  in  the  signal 

A  ciq)ability  of  the  WT  that  has  not  been  discussed  is  noiseAransient  suppression.  Given  a 
signal  that  has  been  ctnrupted  by  some  sort  of  stochastic  process,  it  is  possible  to  deoHnpose  the 
signal  in  a  wavelet  basis  and  then  reconstruct  the  signal  with  the  noise  suppressed.  The  results 
section  concludes  with  an  illustration  of  this  capability.  The  original  signal  f(x)=sin(x)  +  N(0,1)  is 
displayed  in  Figure  9.  The  wavelet  coefficients  for  a  Daubechies  wavelet  N  =  8  is  plotted  in  Fig¬ 
ure  10. 


HGURE  9.  f(x)  =*  sin(x)  +  N(0,1) 
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HOURE10.DAUBECHZES1«^WAVEI'^  OEFHCIENTS 
F(»  SIGNAL  CV?  HCURE !; 

The  aim  is  to  prune  scxne  of  the  coefficients  of  Uie  noisy  signal,  while  the  ultimate  goal  is 
the  removal  of  those  coefficients  produced  by  the  noise  component  of  the  signal  The  hard  thresh¬ 
olding  approach  of  Donoho  and  Johnstone^  is  used  here.  Ln  this  approach,  the  variation  of  the 
coefficients  is  computed  and  a  given  coefficient  d  is  kept  if  it  satisfies 

Idl^s^iogAf  (9) 

where  M  is  the  number  of  wavelet  coefficients.  Thresholding  can  be  applied  to  any  number  of  lev¬ 
els. 

The  pruned  coefficients  from  the  sinusoidal  signal  are  displayed  in  Hgure  1 1.  The  pruning 
operator  was  iqpplied  across  all  the  levels.  As  can  be  gleaned  from  Hgure  1 1,  the  number  of  coef¬ 
ficients  has  decreased  maiicedly. 
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WbwIH  D>ooinpoiltton  CoHleHna 


HGURE 11.  C(»:FiaENTS  AFIER  PRUNING  TECHNIQUES 
ex’ IX^H(»  AND  JOHNSTCX4E 


The  signal  reconstructed  from  the  pruned  coefficients  is  plotted  in  Figure  12.  As  can  be 
seen«  the  pruning  process  has  removed  a  large  ctm^xment  of  the  noise  firom  the  signal.  In  addition, 
the  reconstruction  process  has  faithfully  rqnoduced  tiie  underlying  structure  of  the  data.  This  type 
ci  noise  suppression  technique  or  one  ^  its  variants  would  have  q>plication  to  the  removal  of  bat¬ 
tlefield  transients,  such  as  a^eiy  blasts,  from  the  helict^ter  signids. 


FIGURE  11  SIGNAL  RECONSTRUCTED  FROM 
PRUNED  COEFFICIENTS 
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SYSTEM  DESIGN  CONSIDERATIONS 


Before  concluding,  some  of  the  system  design  and  in^lementation  issues  will  be  consid¬ 
ered.  Some  ideas  will  be  presented  as  to  how  the  WT  can  be  used  as  part  ot  a  signal-processing 
suite,  and  i9>i;nopriate  hardware  for  its  inqilementation  will  be  discussed.  The  use  of  the  WT  for 
feature  extraction  and  nmse  suppression  has  been  demcmstrated  in  the  previous  sections.  The  pro¬ 
posed  classification  system  uti^s  both  of  these  capabilities,  as  illustrtued  in  Hgure  13. 

After  acquisititm,  the  WT  of  the  signal  is  pofOTmed  using  dedicated  hardware,  which  will 
be  discussed  later.  The  signal  is  then  reconstrucu^  without  transients  and  noise  and  passed  to  die 
haimogram  analysis  section.  It  must  be  noted  diat  the  WT  section  functions  not  cxily  to  clean  die 
input  signal  but  also  to  extract  wavelet  based  features  from  the  signal.  These  features  are  then  fed 
to  the  classification  section  for  fusi<»  with  the  harmogram  features.  In  the  classifier  section  of  the 
system,  the  WT  and  harmogram  features  are  fused,  and  discriminant  analysis  is  performed. 

A  good  otample  of  dedicated  WT  hardware  is  the  AWARE  Wavelet  Transform  Processor 
(WTP)  chip.  This  device  is  citable  of  perfrarming  both  the  one-  and  two-dimensional  wavelet 
transform  in  both  the  forward  and  inverse  directions.  This  chip  fully  implements  2  ^larr),  4,  and  6 
coefficient  transfoimations  and  can  inclement  transformations  requiring  a  longer  supported 
wavelet  via  a  cascading  mechanism.^  chip  operates  at  a  clock  speed  of  30  MHz  and  has  a 
reported  execution  speed  of  33  nsec  per  transform^  coefficient  This  speed  is  fast  enough  to  meet 
tlK  requirements  of  most  acoustic  ai^  image  processing  systems. 


FIGURE  13.  SYSTEM  FLOWCHART 
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CONCLUSION 


Some  preUminaiy  results  have  been  provided  that  detail  die  aji^lication  d  the  WT  to 
acoustk  signal  classification.  Whether  die  betovior  oi  die  wavelet  features  is  a  manifestation  of 
die  periodic  component  oi  the  signal  or  an  as-yet  unidentified  transient  conqwnent  will  be  the 
sutiiject  ci  investigations.  Other  future  work  focus  on  the  multiclass  proUem.  the  use  ctf 
nonorthogonal  wavelet  bases,  and  alternate  techniques  for  feature  extraction.  In  addition,  the  tran¬ 
sient  siqipression  aUlity  oi  die  WT  on  field  cdlect^  artillery  blast  data  will  be  evaluated. 
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